Thermodynamic and Tunneling Density of States of the Integer Quantum Hall 

Critical State 
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We examine the long wave length limit of the self-consistent Hartree-Fock approximation irre- 
ducible static density-density response function by evaluating the charge induced by an external 
charge. Our results are consistent with the compressibility sum rule and inconsistent with earlier 
work that did not account for consistency between the exchange-local-field and the disorder po- 
tential. We conclude that the thermodynamic density of states is finite, in spite of the vanishing 
tunneling density of states at the critical energy of the integer quantum Hall transition. 
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A two dimensional(2D) electron gas in the strong mag- 
netic field limit has an zero-temperature insulator-metal- 
insulator transition in the presence of random disor- 
der. When the strength of the magnetic field sweeps 
past a critical value, the off-diagonal (Hall) conductance 
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changes abruptly by %■ while the diagonal (dissipative) 
conductance rises from zero to a finite value of order 
e 2 /2h and then returns to zero. This transition is often 
referred to as the integer quantum Hall transition. The 
filling factor range of the intermediate metallic phase is 
believed to collapse to a conducting critical point in the 
low temperature thermodynamic limit. According to the 
non-interacting theory of the integer quantum Hall effect, 
this behavior occurs because single-electron states are ex- 
tended at an isolated critical energy near the middle of 
each disorder broadened Landau level. The insulator- 
metal-insulator transition takes place when the Fermi 
level of the 2D electron system and one of the discrete 
critical energies cross. 

Our current understanding of the integer quantum 
Hall transition [0 is far from complete, however, since 
Coulomb interaction between electrons are expected to 
provide a relevant perturbation at the non-interacting 
fixed point In addition there are several discrepan- 
cies between the results of the non-interacting theory and 
the experimental findings which may be due to the role 
of interactions. For example, recent experimental work 
has shown that the tunneling density-of-states (DOS) 
vanishes linearly at the Fermi energy || in sharp con- 
trast to the finite DOS of the non-interacting theory. 
To account for this discrepancy Yang and MacDonald 
Q carried out a numerical study in which disorder was 
treated exactly and Coulomb interaction was described 
by a self-consistent Hartree-Fock (HF) approximation. 
They found a linear Coulomb gap at all filling factors 
of the lowest Landau level, even at the critical energy. 
Wang and Xiong J5j investigated the effects of dynamical 
screening in a systematic nonperturbative resummation 



of the most singular diagrams. They found that the DOS 
is linear at the critical energy, as in the Hartree-Fock the- 
ory, exhibiting quantum Coulomb gap behavior. In spite 
of the qualitative DOS change due to interactions, how- 
ever, Yang et al. |J found that the value of the localiza- 
tion length exponent vi oc is identical to that of the non 
interacting theory ||7|-|i"l|| . 

A second discrepancy between non-interacting electron 
theory and experimental findings is the value of the dy- 
namical scaling exponent, z. For non- interacting elec- 
trons, z equals the space dimension d = 2. This re- 
mains true for short-ranged interactions since they are 
irrelevant at the non- interacting fixed point [§,[l2| . How- 
ever, experimental data are consistent with z Rj 1 [fl3|| . 
Lee and Wang j| conjectured that a change in z with 
no change in vi oc is due to the non-critical linear sup- 
pression of the single-particle DOS induced by Coulomb 
interactions. [§,1). Recently, Huckestein and Backhaus 
fLU attempted to substantiate this conjecture by evalu- 
ating the density-density response function in the quan- 
tum Hall regime, including interaction effects within a 
time dependent Hartree-Fock approximation (TDHFA). 
Their analysis of dynamic scaling gives z — 1 and ap- 
pears to be consistent with the result of a naive scaling 
argument incorporating the linearly vanishing DOS: the 
frequency-dependent length scale L w — uj^ 1 ^ = , 1 

with density-of-states po cx u> implies z = 1. However, it 
is not clear that the po that enters L w is the tunneling 
DOS. In fact, it has been argued jl5| that the relevant 
DOS here is the thermodynamics DOS or the compress- 
ibility dn/dp, which should be smooth and finite for a 
disordered system on general grounds. Huckestein and 
Backhaus, however, found that the thermodynamic DOS 
or the compressibility is zero at the transition. Since the 
compressibility is proportional to the inverse screening 
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2ne ^j, their result would imply an in- 



length, q sc 

finite screening length. The screening properties of the 
conducting critical state would therefore be highly uncon- 



ventional. This result is plausible since the single particle 
states at a quantum Hall transition have the special prop- 
erty that only states at the critical energy are delocalized 
while all other states are localized. However, in order to 
make their calculation numerically manageable, Huck- 
estein and Backhaus [H were forced to neglect corre- 
lations between the Hartree-Fock theory self-energy and 
the random potentials. Closely related approximations 
have long been common in a number of circumstances, 
for example in addressing the properties of disordered su- 
perconductors where it has recently been demonstrated 
that they are unreliable. ]l6[ |. 

In order to evaluate the Hartree-Fock irreducible 
density-density response function directly, it is neces- 
sary to solve a complicated integral equation. Huckestein 
and Backhaus circumvented this difficulty by disorder- 
averaging the quasi-particle response function and the 
self-energy separately. We circumvent this difficulty 
without making approximations, by evaluating only the 
total charge response Q to the external potential from 
a point charge, solving the Hartree-Fock equations for 
the self-energy that is consistent Jjj]] with each disorder 
realization. In linear response, Q = lim^o x(9)27re 2 /eg 
where x(q) is the full static response function, related 
to the irreducible response by x(q) = X %rr (<z)/[l + 
(2we 2 / eq)x irr (q)}, and e is the dielectric constant. It fol- 
lows that 
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q ^o 2i:e 2 x lrr {q) 
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In the approximation studied by Huckestein and Back- 
haus Jl4| , it was determined numerically (see their Fig. 3 
and discussions) that XhbG?) = s ^? + 0{£/L), where s ~ 
0.2 is a numerical constant and I is the magnetic length. 
This would imply that the second term on the right hand 
side of Eq. (Q) has a finite value and Qhb — 0.56. On 
the other hand, it is normally the case for metals that 
the screening wave- vector q sc = 2ire 2 lim 9 _>o X trr {l)i is 
related to the thermodynamic density-of-states according 
to q sc — 2ne 2 dn/dfj, ]lq| , a result often referred to as the 
compressibility sum rule, which would imply that Q = 1. 
(It is clear from experiment that dn/dfi must be finite at 
the integer quantum Hall transition.) Our numerical re- 
sults support the applicability of the compressibility sum 
rule and are clearly inconsistent with the value of Qhb 
that violates this sum rule by about 50%. When this 
rule is satisfied Q = 1, the external charge is perfectly 
screened, and the screening length q sc is finite. Our re- 
sults show that despite the linearly vanishing Coulomb 
gap in the tunneling DOS, x irr (<? = 0) is finite and im- 
ply that transport is indeed governed by diffusion, as for 
non-interacting electrons. 

We perform our calculations in the Landau gauge 
[A = (0, Bx, 0)] and apply quasi-periodic boundary con- 
ditions to the HF single-particle orbitals inside a square 
with area A — a 2 . The basis states used to represent the 



HF Hamiltonian are related to elliptic theta functions 
and can be labeled by a set of guiding centers \Xj >, 
j = 1,...,JV^ inside the fundamental cell of the finite 
system is the total flux quantum passing through 
the area A). The wave functions of the basis states are 
< r\Xj >= IX-oo faj+kaif) where 



<l>x{r) = 
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(2) 



The HF eigenstate is given by a linear combination of the 
basis states. \a >— J2i \Xi >< Xi\a >. The coefficients 
< Xi\a > satisfy a set of matrix equations for the HF 
Hamiltonian, Hhf = Vhf + Vj + Vi on , which is a sum 
of the HF potential in the LLL, the impurity potential, 
and the ionic potential of the inserted test charge: 

Y}< X^Xi > + < x y \v wn \x 3 > 

i 

+ < Xj\V H F\Xi >] < Xi\a >= e a < Xj\a > . (3) 

The position of the ion is R = a(i, 5). The matrix ele- 
ments of the ion potential is 

1 ^ 2vre 2 



< Xji\Vi on \Xj >= 



—T — 

A eq 
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x6' j/d+my exp(i fe (X, + q y £ 2 /2) - q 2 £ 2 /4), (4) 

where q = 2 f{m x + sN^, rriy + tN^) with m x>y = 1, 
and s,t are integers. We have used a model disorder 
potential consisting of Nj delta-function scatters with 
strength uniformly distributed between —A and A and 
scattering centers uniformly distributed inside the fun- 
damental cell of the finite-size system. Hartree and Fock 
potentials can be expressed as a function of the electron 
density 



< X s \VHF\Xi > =J2 A (q) ex P [ 



iq^Xi + Xj). 



x5'(j,i + m y )U H F(q)- 



(5) 



Here the sum over q is over the discrete set of wave vectors 
consistent with the boundary conditions and S'(n,n') is 
1 if n' = n(modA^0) and otherwise. Uhf{<i) is propor- 
tional to e 2 /e£ and includes both Coulomb and exchange 
interactions. The quantity A(q) is proportional to the 
Fourier component of the charge density and is calcu- 
lated self-consistently from the eigenvectors of the HF 
Hamiltonian by 



xexpl 



v j=i y=i 
MXi + Xj) 



N 



]J2<Xf\a><a\X J >, (6) 
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where N is the number of electrons in the system. For 
each disorder realization the matrix equation is solved for 
several values of 7 = (e 2 /e£)/T where T = (\ 2 N I /£ 2 A) 1 / 2 
is the characteristic disorder potential energy scale and 
Ni is the number of impurities. For this model the self- 
consistent Born approximation (SCBA) DOS is nonzero 
for |e| < (2/37r)( 1 / 2 )r - 0.46r. 

For a system of N electrons in the fc-th realization of 
the random disorder potential, we write the electron den- 
sity as p%{r) = £ a =i,...,iv |*£ \r)\ 2 where P = 0,1, 
depending on whether the ion is absent or present. Let 
Po,n be the density of the uniform background charge. 
When Pq N (r) is disorder-averaged over No number of re- 
alizations we find po,N(r) = jj^ J2k Pojfi^) = Po,N- The 
charged density induced by the ion may be evaluated 
by co mputing P/,j y+i(r) - p 0:N {r) = pi,N+l(r) ~ Po,N- 
Since po,jv(^0 = Po,N+i(r) — 4 the induced density 
Ap(r) = p I<N+1 (r) - p ,N+i{r) + jj. 

Figure 1 displays the charge induced when an ion is 
placed at the center of the square at LLL filling fraction 
v = 1/2 and for the parameters 7 = 0.6, No — 4 and 
N,p — 200. The inset in Fig. 1 shows how the total in- 
duced charge changes as the distance from the ion is var- 
ied. These results are consistent with Q = 1 and a finite 
screening length A sc ~ hi. They are clearly inconsistent 
with Huckestein and Backhaus's result Qhf — 0.56. Us- 
ing the DOS of the SCBA, we obtain a shorter Thomas- 
Fermi screening length, X TF = 2ir/(je) ~ 11. Fig- 
ure 2 shows the disorder averaged Ap(x,yi) for a much 
broader Landau level with 7 = 0.1, v = 1/2, No = 10, 
and N$ = 162. In this case we find that the screen- 
ing length X sc — ll£ while Xtf — 61 Our results for 
the 7-dependence of the screening length is in qualitative 
agreement with that implied by the DOS behavior in the 
SCBA. For a more quantitative comparison for the dis- 
order dependence of the screening length, a Hartree non- 
linear screening treatment, in which disorder is treated 
in the SCBA p9| , may be required instead of the simple 
SCBA. 

We have also investigated the screening properties in 
the insulating regime at v = 1/7. Figure 3 displays the 
induced densities for 7 = 0.1, No = 10, and = 175. 
We observe from the plot that the induced density does 
not decay rapidly and that the fluctuations are spread 
over a wide range. In this case it is difficult to define a 
screening length. We believe this reflects the incomplete 
screening property pC| of an insulator. Indeed from the 
inset we see that the screening length appears to diverge 
with the finite system size. However, in the presence of 
a weaker random disorder with 7 = 0.6, where the ten- 
dency towards crystallization is strong, we find a different 
result. Contour plots of Po°\r) and p^(f) are shown in 
Fig. 4. Plot (a) is without the ion while plot (b) is with 
the ion. Comparing these two plots we can infer that the 



ion captures an electron while the location of the other 
electrons are unchanged. A similar abrupt reconfigura- 
tion of the ground state was found and studied in the 
pinning of Wigner crystals [§l|. 

We have shown here that despite the linearly vanish- 
ing Coulomb gap in the tunneling DOS, the thermody- 
namic DOS of the integer quantum Hall critical state 
remains finite, and transport is indeed governed by dif- 
fusion as for non-interacting electrons. It still remains 
an open question whether an improved implementation 
of the TDHFA, that takes into account self-consistently 
self-energy and vertex corrections, can give rise to the ob- 
served dynamic critical exponent for the integer quantum 
Hall transition. 
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FIG. 1. Induced charge distribution f or the c riti- 
cal state. Each curve represents one of Ap(x, yC) for 
i = 10, 11, 19, 20, where Vi = ^T 1 . These curves are 
shifted vertically for the sake of clear presentation. The 
parameters are v — 1/2, 7 = 0.6, No = 4, and iV^ = 200. 
The unit of length is I. Inset: Integrated induced charge 
as a function of the distance to the position of the ion. 
The dashed line corresponds to Qhb — 0.56. 
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FIG. 2. Same as in Fig.l but for a more disorder 
broadened LLL with 7 = 0.1, No — 10, and = 162. 
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FIG. 3. Same as in Fig.l but for an insulating state 
with v = 1/7, 7 = 0.1, N D = 10, and N^ = 175. 
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FIG. 4. Contour plots of the charge density in the pres- 
ence of a strong the tendency towards crystallization for 
7 = 0.6. at v = 1/7 with N D = 1 and = 175. Plot (a) 
is without the ion while plot (b) is with the ion. Com- 
paring these two plots we can infer that the ion captures 
an electron while the location of the other electrons are 
unchanged. The unit of length is I. 
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